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ABSTRACT  OF  OBJECTIVES  AND  ACCOMPLISHMENTS 

The  goal  of  this  research  program  was  to  consider  the  electronic  structure  aspects  of  radiative 
and  radiationless  decay  processes  related  to  the  stability  and  formation  of  high  energy  density 
materials.  In  order  to  accomplish  this  goal  a  unique  system  of  electronic  structure  algorithms,  the 
BROOKLYN  programs,  has  been  developed.  These  programs  provide  advanced  capabilities  for  the 
study  of  the  electronic  structure  aspects  of  spin- forbidden  and  electronically  nonadiabatic  processes 
The  methodology  which  we  developed  is  based  exclusively  on  large  scale  configuration  state  function 
expansions  (105— 106  terms).  These  methods,  which  we  believe  define  the  state  of  the  art  in  these 
areas,  have  permitted  us  to  make  significant  contributions  to  the  understanding  of  radiative  and 
radiationless  decay  processes.  Problems  of  particular  relevance  to  the  high  energy  density  materials 
program  include:  (i)  a  study  of  the  possible  stability  of  the  energetic  maximum  ionicity  state  (HeH)+- 
H-;  (ii)  a  study  of  the  spin-forbidden  decay  mechanism  of  the  model  azide  system  N3H(1A')  — > 
N2(1Lg)+  NH(X3E‘)  and  finally  a  study  of  the  decay  pathways  for  the  metastable  excited  state  of  the 
helium  atom,  He(23S),  resulting  from  collisions  with  other  (ground  state)  helium  atoms. 

Finally  a  new  phase  of  program  development  has  been  initiated  focusing  on  the  efficient 
determination  of  regions  of  allowed  and  avoided  crossings  of  potential  energy  surfaces  using  analytic 
gradient  techniques.  This  program  development,  which  will  be  completed  during  our  present  AFOSR 
sponsored  research  program,  will  provide  powerful  new  tools  for  the  study  of  spin- forbidden  and 
electronically  nonadiabatic  processes. 
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Stability  of  Potential  High  Energy  Density  Materials 

Final  Technical  Report 
(1  July  1986-30  October  1989) 


SUMMARY 


The  goal  of  this  research  program  was  to  consider  the  electronic  structure  aspects  of 
radiative  and  radiationless  decay  processes  related  to  the  stability  and  formation  of  high  energy 
density  materials.  In  order  to  accomplish  this  goal  a  unique  system  of  electronic  structure 
algorithms,  the  BROOKLYN  programs,  has  been  developed.  These  programs  provide  advanced 
capabilities  for  the  study  of  the  electronic  structure  aspects  of  spin-forbidden  and  electronically 
nonadiabatic  processes  The  methodology  which  we  developed  is  based  exclusively  on  large  scale 
configuration  state  function  expansions  (10^-10^  terms).  These  methods,  which  we  believe  define 
the  state  of  the  art  in  these  areas,  have  permitted  us  to  make  significant  contributions  to  the 
understanding  of  radiative  and  radiationless  decay  processes.  Problems  of  particular  relevance  to 
the  high  energy  density  materials  program  include:  (i)  a  study  of  the  possible  stability  of  the 
energetic  maximum  ionicity  state  (HeH)+-H*;  (ii)  a  study  of  the  spin-forbidden  decay  mechanism 
of  the  model  azide  system  N3H(1A')  -»  N2(1£^)+  NH(X3Z_)  and  finally  a  study  of  the  decay 

o 

pathways  for  the  metastable  excited  state  of  the  helium  atom,  He(23S),  resulting  from  collisions 
with  other  (ground  state)  helium  atoms. 

Finally  a  new  phase  of  program  development  has  been  initiated  focusing  on  the  efficient 
determination  of  regions  of  allowed  and  avoided  crossings  of  potential  energy  surfaces  using 
analytic  gradient  techniques.  This  program  development,  which  will  be  completed  during  our 
present  AFOSR  sponsored  research  program,  will  provide  powerful  new  tools  for  the  study  of 
spin-forbidden  and  electronically  nonadiabatic  processes. 
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TECHNICAL  REPORT 


One  of  the  principal  aspects  of  our  research  program  has  been  the  development  of  new 
theoretical  tools  for  treating  spin-forbidden  and  electronically  nonadiabatic  processes.  These  tools 
have  then  been  used  to  study  the  stability  and/or  decay  mechanisms  of  potential  high  energy 
density  materials(HEDM's).  We  believe  that  as  a  result  of  computational  limitations  the  electronic 
structure  aspects  of  these  important  classes  of  chemical  problems  had  received  inadequate 
treatment.  The  methodology  which  we  developed  is  based  exclusively  on  large  scale  configuration 
state  function  expansions  ( 1 05—  1 06  terms).  These  methods,  which  we  believe  define  the  state  of 
the  an  in  these  areas,  have  permitted  us  to  make  significant  contributions  to  the  understanding  of 
radiative  and  radiationless  decay  processes.  Below  we  outline  the  methods  we  have  developed  and 
their  applications  in  the  HEDM  program 
I:  SPIN-FORBIDDEN  CHEMISTRY 

Our  studies  in  the  area  of  spin-forbidden  chemistry  employ  the  Breit-Pauli  approximation1 
and  consider  those  interactions  which  lift  the  2S+1  degeneracy  associated  with  the  total  electron 
spin  quantum  number,  S.  These  interactions  are  the  full  microscopic  spin-orbit  interaction  A80, 

h“(i)  s,  +  X  h‘"(j,;)  (s,+2sJ) 

i=l  i,j= 1  1  •  1 


and  the  dipolar  spin- spin  interaction  A® 


S  s  .  - hss(  •  s. 

i,j=\ 


where  h^fi),  hsoo(i,j),  and  h^fij)  are  the  microscopic  spin-orbit,  spin-other-orbit,  and  dipolar 
spin-spin  interactions1  and  we  write  ll Bp  =  $ 80  +  M. 

(A)  Methodology 
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Our  treatment  of  these  interactions  is  unique  in  two  ways.  We  use  the  symbolic  matrix 
element  methods  originally  introduced  in  the  treatment  of  the  nonrelativistic  Bom  Oppenheimer 
hamiltonian  A0  by  Liu  and  Yoshimine2  to  evaluate  the  matrix  elements  of  and  Ass4  in  the  N- 
electron  configuration  state  function  (CSF)  space.  This  enables  the  use  of  large  scale  CSF  spaces 
( 1 05—  1 06  terms)  to  characterize  these  interactions.  These  spaces  are  one  to  two  orders  of 
magnitude  larger  than  those  accessible  using  previously  existing  methods1'6  and  have,  as 
discussed  below,  improved  significantly  the  reliability  of  theoretical  work  in  this  area. 

The  second  unique  aspect  of  our  approach  involves  the  manner  in  which  the  relativistic 
wavefunction  is  obtained  from  the  zeroth  order  nonrelativistic  [configuration  interaction  (Cl)] 
wavefunction.5  Within  the  context  of  the  Breit-Pauli  approximation  the  relativistic 
wavefunction.'Ffl),  is  obtained  from  the  zeroth  order  nonrelativistic  wavefunction,  VP°(I),  using 
perturbation  theory.  The  commonly  used  approach  had  been6  to  determine  the  perturbation 
contribution  to  the  relativistic  wavefunction 


^(1)  =  'FO(I)  +  ¥1(1) 


1.3 


from  first  order  perturbation  theory  using  a  spectral  resolution 


^\l)  =x 

K 


(y°u) 


H 


BP  v 


E°U)-  E°U) 


1.4 


where  ^(K)  satisfies  the  nonrelativistic  Schroedinger  equation 

[AO-EO(D]vFO(I)  =  0 


1.5 
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This  procedure  is  unsatisfactory  owing  to  the  need  to  drastically  truncate  the  sum  in  eq.  1.4. 

Instead  we  have  introduced5  a  technique  which  allows  VF1(I)  to  be  determined  exactly,  within  a 
given  CSF  space,  as  the  solution  of 

[ri0-E0(I)]4/l(I)=-ABP4'0(I)  .  1.6 

The  use  of  eq.  1.6  in  lieu  of  eq.  1.4  has  improved  considerably  our  ability  to  predict  radiative  rates 
for  spin-forbidden  dipole-allowed  processes. 

Eq.  1.6  is  only  satisfactory  when  the  electronic  states  in  question  are  well  separated.  When 
one  or  more  of  the  electronic  states  in  eq.  1.4  is  near  to,  quasidegenerate  with,  ^(I)  a 
generalization  of  eq.  1.6  is  required.  As  part  of  this  research  program  we  have  shown7  how  such 
a  generalization  can  be  derived  from  partitioning  theory.^  The  description  which  emerges  from  this 
formalism  expresses  the  electronic  wavefunctions  as  linear  combinations  of  dressed  diabatic 
wavefunctions  ['Fd(I)] 


4'd(/)  =  +  L7 

where  the  dressing  term,  VP<^ d  (/)  represents  the  perturbation  of  the  (diabatic)  quasidegenerate  state 
VF°(I)  by  states  outside  the  quasidegenerate  space  and  satisfies 


[H°  -  £°(/)]'Pq<,(/)  =  -QHBPx¥°{I) 


1.8 


where  0  projects  onto  the  orthogonal  complement  of  the  quasidegenerate  space.  The 
quasidegenerate  dressed  diabatic  states  xPd(I)  which  are  the  analogues  of  the  'F(I)  in  eq.  1.3  are 
mixed  by  the  Breit-Pauli  interaction  ABP.  This  mixing  is  given  by  a  secular  equation  with  the 
dimension  of  the  quasidegenerate  space  which  has  the  form: 
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1.9 


^1,  1  E  ^1,2 

1  H2,2~E 


4^(1) 


where  H,  y  =  ^vF‘i(/)|//  °+  H  BPJvF‘f(y)y  .  Since  the  coupling  between  the  states  4/d(I)  originates 
from  a  potential  term  ltBP  rather  than  the  derivative  coupling  terms  to  be  discussed  subsequently 
the  4/d(I)  are  referred  to  as  diabatic  states. 

Applications  of  the  above  described  methodology  in  the  area  of  spin-forbidden  radiative 
decay  are  discussed  below. 

(B)  Applications 

(i)  Radiative  Decay  in  Nitrogen  Halides,  NF  and  NCI.  Resolution  of  Significant  Experimental 
Discrepancies 

The  lifetime  with  respect  to  spin-forbidden  radiative  decay  is  an  important  attribute  when 
assessing  the  utility  of  a  metastable  species  as  an  energy  transfer  agent  in  a  chemical  laser  system. 
The  b1^  and  a!A  states  of  the  nitrogen  halides  NF  and  NCI  are  such  metastable  states  as  in  each 
case  the  molecule  has  an  X3Z'  ground  state.  These  systems  have  therefore  been  suggested  as 
potential  energy  sources.  This  possibility  motivated  our  theoretical  studies  of  the  spin-forbidden 
dipole-allowed  radiative  decay  processes  (b1£+,  a1  A)  ->  X3Z_  in  both  NF  and  NCI. 

Our  study  of  the  b!Z+  — >  X3Z‘  transition  in  NF  revealed  a  key  misconception  concerning 
the  nature  of  this  transition.  It  had  been  previously  reported9  that  pii/pj.  ~  1/7  based  on 
Watson's10  theoretical  model  of  the  experimental  results  of  Douglas  and  Jones  (DJ).11  Here 


1.10a 
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1.10b 


^=('K*1ztJ|£iK*3Z7)) 


where  is  the  electronic  dipole  moment  operator  and  'F(I)  (eq.  1.3)  was  obtained  from  eqs.  1.5- 
1 .6.  Although  this  result  had  been  confirmed  by  SCF  level  calculations9  it  is  in  fact  in  error! 

Using  the  basic  data  of  DJ  and  Watson's  theoretical  model  we  observed5  that  the  correct  ratio 
should  be  ~  7,  that  is  the  inverse  ratio  had  been  reported!  As  a  result  of  the  confusion  over 
the  original  experiment  of  DJ,  Vervloet  and  Watson12  reinvestigated  the  0-0  band  of  the  b1^  — > 
X3Z~  emission  using  a  Fourier  transform  spectrometer.  Their  result  fJ.il/Ma  =  3.8  is  in  excellent 
agreement  with  our  best  theoretical  estimate13  lUi/jii  ~  3.87. 

For  the  a1  A  state  of  NCI  the  situation  was  much  less  subtle.  In  this  case  conflicting  lifetime 
data  for  the  a1  A  — » X3£-  transition  had  been  reported.14-15  Coombe  and  van  Benthem143  had 
reported  a  lifetime  of  205  nsec  based  on  gas  phase  measurements.  However  shortly  thereafter 
Becker  and  coworkers15  reported  a  lifetime  based  on  matrix  isolation  experiments  of  1440  nsec 
(which  extrapolates  to  a  gas  phase  value  of  2658  nsec).  Our  calculations  (reported  as  a 
communication)16  gave  an  estimated  lifetime  for  the  a1  A  -»  X3L~  transition  of  2400  nsec 
confirming  the  results  of  Becker  et  al. 

(ii)  Spin-Forbidden  Radiative  Decay  Resulting  from  Coupling  to  Bound  States  Embedded  in  a 
Continuum.  Application  to  the  a1  A  — ¥  X  32r  Transition  in  CH~ 

The  previously  noted  work  on  the  nitrogen  halides  not  withstanding  perhaps  the  most 
convincing  demonstration  of  the  power  of  the  methodology  we  have  introduced  comes  from  our 
treatment17  of  the  aJA  — >  X32r  transition  in  CH_  a  transition  which  had  been  studied 
experimentally  by  Lee's  group  at  Berkeley.18  This  transition  borrows  intensity  from  3n,  *11  states 
embedded  in  the  ionization  continuum  of  CHf^ri)  +  e~  Consequently  traditional  techniques  for 
handling  this  problem  based  on  the  eigenstate  expansion  of  ^(I),  eq.  1.4,  are  computationally 
intractable  because  of  the  necessity  of  obtaining  a  large  number  eigenstates,  'PO(K)  of  A0. 
Consequently  Lee  and  coworkers  note  in  their  work  that  a  theoretical  estimate  of  the  lifetime  of  the 
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a1  A  — »  X^Ir  f  ansition  in  CH_  "would  be  interesting  but  difficult  given  the  presence  of  the 
ionization  continuum."  However  our  approach  overcomes  this  difficulty  since  by  obtaining  ^'-(I) 
directly  from  eq.  1.6  the  correct  contribution  from  each  eigenstate  is  obtained. 

A  related  problem  concerns  the  determination  of  a  molecular  orbital  space  appropriate  for 
the  characterization  of  a  function  which  contains  contributions  from  highly  excited  states. 
We  have  developed  a  method  to  address  this  question  based  on  the  iterative  natural  orbital  (INO)19 
procedure. 

Our  treatment  of  the  spin-forbidden  radiative  transition  a!A  — »  X^L-  in  CH~  has  been 
reported.17  Using  the  INO  procedure  it  was  found  that  the  optimum  orbital  space  for  describing 
^(I)  includes  a  molecular  orbital  with  character  intermediate  between  the  compact  valence  orbital 
and  the  diffuse  orbital  obtained  from  two  alternative  MCSCF  procedures.  This  orbital  space,  the 
INO  orbital  space,  was  obtained  from  these  two  distinctly  different  starting  points  providing 
convincing  evidence  of  the  utility  of  this  approach  for  obtaining  molecular  orbitals  for  use  in  the 
description  of  T^I).  Using  the  INO  orbital  set  a  lifetime  (x)  for  the  ground  vibrational  level  of  the 
aJA  state  of  6.14(±1.2)s  was  obtained.  This  value  is  in  excellent  agreement  with  the  experimental 
value  t  =  5.9(+0.8,  -0.6)s  reported  by  Okumura  et  al.1**  Note  that  agreement  with  experiment 
would  have  been  less  satisfactory  (x  =  4.41s)  if  the  INO  optimized  orbitals  had  not  been  used. 

(iii)  Spin-Forbidden  Radiative  Decay  Involving  Quasi-Degenerate  States.  Application  to  the 
B1£+  — »  a3n  Transition  in  MgO. 

As  noted  above  when  treating  spin-forbidden  radiative  processes  involving  spin-allowed 
surface  crossings  or  closely  spaced  PES's,  that  is  spin-forbidden  radiative  decay  involving 
quasidegenerate  states,  the  straightforward  first  order  perturbation  treatment  given  by  eqs.  1.4  or 
1 .6  is  inappropriate.  Instead  the  method  outlined  in  eqs.  1.7- 1.9  is  required.  This  approach 
recognizes  the  special  status  of  the  quasidegenerate  states  while  retaining  the  advantages  inherent  in 
the  use  of  eq.  1 .6  (rather  than  eq.  1.4)  for  treating  the  perturbation  resulting  from  the  remaining 
states. 
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Quaside generate  perturbation  theory  was  used  to  study7  the  spin-forbidden  dipole-allowed 
B1!4-  — »  a3n  transition  in  MgO  which  had  been  investigated  experimentally  by  Ip  in  R.  W.  Field's 
group  at  MIT.20  In  order  to  characterize  this  transition  five  (dressed)  diabatic  wavefunctions  were 
determined 

'F(/4^)  =  'F0(,4Ini) 

¥(a3rj)  =  +  4'g(1ni;a3ni) 

4'(x,s;.)  =  >p“(x,2;;,)+y(X;xlz;,)+^(In(,.;x,x;.) 

'•’(a’H.) = *  Vn.) + 

where  the  subscript  Q  indicates  that  the  first  order  perturbation  contribution  (the  dressing)  is 
determined  in  a  subspace  orthogonal  to  the  lowest  state  of  the  indicated  symmetry  which  in  the 
present  circumstance  is  quasidegenerate  with  the  zeroth  order  state  being  perturbed.  The  dressing 
terms  OF1)  give  rise  to  an  ’intrinsic'  spin-forbidden  transition  moment  for  the  B!L+  — >  a3n 
transition  as  they  result  in  wavefunctions  of  mixed  singlet-triplet  character.  The  geometry 
dependent  mixing  of  the  pair  of  diabatic  states  in  eq.  1.12  and  (separately)  the  pair  in  eq.  1.13  is 
determined  from  a  secular  equation  with  hamiltonian  1^°  +  so  (eq.  1 .9).  This  mixing  gives  rise 
to  an  additional  contribution  to  the  B1I+  — » a3n  transition  moments. 

Using  this  approach  the  parallel  wKB1^.,  v  =  0;  a  3rio+,  v  =  n)  and  perpendicular 
v  =  0;  a  3rii,  v  =  n)  components  of  the  B1L+  a3n  transition  moment  were 
determined.  The  n  =  0, 1  transitions  had  been  observed  experimentally.20  The  results  for  these 
transitions  were  found  to  be  in  good  accord  with  the  available  experimental  data.  For  example  the 
calculated  ratio  F  of  the  radiative  rate  for  the  B1!4-  — >  a3n(0,0)  transition  and  the  B1^4"  -> 
a3n(0,0)  transition  is  F  =  5.37  x  lCh3  which  compares  well  with  the  estimate  inferred  from  the 
work  of  Ikeda  et  al.21,  F  =  5.6  x  10~3. 

(iv)  Lifetimes  of  Metastable  (Potential  Energetic)  Species:  The  Quenching  of  Helium  2^S 


1.11 

1.12a 

1.12b 

1.13a 

1.13b 


This  study22  was  concerned  with  the  lifetime  of  the  metastable  energetic  species  He(23S)  in 
liquid  helium  considering  the  possible  role  of  nonadiabatic  effects  in  the  collisional  quenching  of 
this  speciesAvhich  energetic  by  1  leV  with  respect  to  the  ground  state  atom.  The  electronic 
structure  aspects  of  a  nonadiabatic-radiative  decay  mechanism  were  analyzed.  The  basic  issue  was 
the  lifetime  of  the  a3E^"  state  of  He2  which  correlates  with  He(23S)  and  ground  state  helium.  In  the 
proposed  mechanism  the  a3X^  state  is  coupled  by  relativistic,  rotational  and  radiative  interactions 
to  the  A1  X*  state  which  serves  as  a  gateway  to  the  XJX+  (electronically  quenched)  state  of  He2 

5 

through  the  spin-allowed  dipole-allowed  bound-free  transition  A!X^  — »  X!Xg. 

State  averaged  MCSCF/second  order  Cl  wavefunctions  for  the  ground  X'X*  state,  and  the 

o 

excited,  a3X^,  AlX„,  b3ng,  Blng,  c3X*,  and  C^X*  states  (referred  to  here  as  the  primary  space) 

o  o 

of  He2  were  determined.  Using  these  wavefunctions  all  interstate  matrix  elements  of  the  form 
(T'OfJ)  I  Ol  VT°(I))  were  determined  for  (i)  O  =  HBP  =  Hso  +  Hss  where  Hso  and  Hss  are 
respectively  the  spin-orbit  and  dipolar  spin-spin  interactions  in  the  Breit-Pauli  approximation,  (ii) 

AAA  AAA 

0=  IS,  where  15  is  the  total  electronic  orbital  angular  momentum  operator  and  (iii)  0=  p.  where  p. 
is  the  dipole  moment  operator. 

In  the  nonrotating  molecule  these  interactions  give  rise  to  the  spin-forbidden  dipole-allowed 
radiative  transitions  (b3ng,  c3X*)  — >  AJX*.  However  a  complete  description  of  these  radiative 
decay  processes  requires  consideration  of  interactions  originating  outside  the  primary  space.  Thus 
in  this  work  the  spin-forbidden,  dipole-allowed  perpendicular,  p^(J,  A^+q+X  J  =  c3Xgj,  b3ngi 
and  parallel,  pn(b3Ilgo+»  AJX^+),  transition  moments  were  determined  using  quasidegenerate 
perturbation  theory. 

The  computed  potential  energy  curves,  coupling  matrix  elements  and  dipole  moments 
permit  a  fully  quantum  mechanical  analysis  of  the  nonadiabatic-radiative  quenching  mechanism. 
However  a  preliminary  phenomenological  analysis  of  aspects  of  this  process  based  on 
quasidegenerate  perturbation  theory  indicated  that  this  mechanism  would  not  be  competitive  with 
the  direct  spin-forbidden  radiative  decay  a3!*— ►X1!*  studied  previously  by  Chabalowski  et.  al.23 

u  o 

(v)  Spin-forbidden  Radiationless  Decay:  N3H(X1A')^N2(X1I^)+NH(X3E-) 

© 
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Recently  there  has  been  considerable  interest  in  the  decomposition  pathways  of  the 
homologous  azides  N3-X  for  X=H 24-26,  p  21  ^  Q  28-29  and  gr  29-30^  that  is 

N3X(X1A*)  -)•  NXCX^Z-a^)  +  N2(1Ig)  1.14 

as  well  as  in  the  isovalent  systems  FNCO  and  NNOO31.  Interest  in  these  systems  has  been 
motivated  in  part  by  the  potential  of  N3F  (fluorine  azide)  as  a  clean  source  of  NF^A).  The 
approximately  1.5eV  of  energy  in  NF^A)32  would  then  be  available  for  use  in  a  collisional 
energy  transfer  chemical  laser  system276  analogous  to  the  oxygen-iodine  laser  system  system.33 
Several  recent  experimental  and  theoretical  studies34-3*’  have  focused  on  reaction  1  for  X=H,D  and 
contributed  significantly  to  a  detailed  understanding  of  this  reaction.  The  production  of  NH(X3Z_) 
in  reaction  1  is  clearly  spin-forbidden.  Alexander,  Werner  and  Dagdigian  (AWD),26  in  a  recent 
theoretical  study  of  the  N3H(hydrazoic  acid)  decomposition,  noted  three  factors  as  beirg  essential 
to  the  understanding  of  the  electronic  structure  aspects  of  this  process:  (i)  the  location  of  the 
minimum  energy  singlet-triplet  crossing  which  in  this  system  represents  the  barrier  to  the  spin- 
forbidden  reaction,  (ii)  the  gradients  of  the  singlet  and  triplet  potential  energy  surfaces  at  this 
crossing  and  (iii)  the  spin-orbit  induced  singlet-triplet  coupling  in  the  vicinity  of  the  minimum 
energy  crossing.  Previous  ab  initio  studies,  including  the  study  of  AWD  which  motivated  the 
present  work,  had  provided  considerable  insight  into  the  location  and  topology  of  the  minimum 
energy  crossing.  However  there  had  been  no  treatment  of  the  spin-orbit  induced  coupling  of  the 
lowest  singlet  and  triplet  states  in  this  key  region  using  extended  Cl  wavefunctions.  To  address 
this  deficiency  we  undertook  a  study  of  the  electronic  structure  aspects  of  the  spin-forbidden 
decomposition  reaction  N3H(X  1A')  — >  NH(X3!-)  +  N2(12^)  37.  Using  a  flexible  basis  of  better 
than  double  zeta-polarization  quality  and  Cl  expansions  of  approximately  200,000  terms  the  spin- 
orbit  (H80)  couplings  between  the  lowest  singlet  state  4/ia,(1A')  s  T'fUA'fO)]  and  the  components 
of  the  lowest  triplet  state  4'ia-(3A")  2  i  'F[13A"(0)],  vF2a<3A")  =  i  {^[13A”(D]  -  'f'[l3A"(-l)]}/V2 
were  determined  in  the  asymptotic  region  corresponding  to  N2+NH,  at  the  experimental36  (planar) 
equilibrium  geometry,  and  in  the  vicinity  of  the  approximate  minimum  energy  crossing.  At  the 
approximate  minimum  energy  crossing  we  finct  h^°  2  OFia'OA')  I  H80 1  vI/ia-(3A"))  =  39.5  cm-1  » 
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hy°  =  OPia'^A')  I  Hs0  I  H/2a'(3A")>  *  0.456  cm*1.  The  negligibility  of  h^°  is  attributed  to  the 
geometry  of  the  N3H  moiety  at  the  approximate  minimum  energy  crossing  and  the  absence  of 
NH^Il)  character  in  the  1  !A'  state.  The  efficacy  of  h^°  in  promoting  the  spin-forbidden  reaction 
1  was  estimated  using  a  Landau-Zener  model.38  It  was  concluded  that  the  spin-forbidden  process 
would  be  significant  only  near  threshold.  This  conclusion  was  attributable  to  both  the  small 
magnitude  of  h*°  and  the  comparatively  large  gradient  of  AE  with  respect  to  the  approximate 
reaction  coordinate  R(N2-N3). 


II:  ELECTRONIC  NONADIABATICITY  INDUCED  BY  DERIVATIVE  COUPLINGS 
The  decay  of  an  electronically  excited  energetic  species  may  involve  a  radiationless 
pathway.  This  pathway  would  in  turn  involve  nuclear  motion  on  more  than  one  adiabatic  Bom- 
Oppenheimer  potential  energy  surface  (PES)  with  the  coupling  between  the  surfaces  arising  from 
derivative  couplings  of  the  form 


gau,  /,  R)  =  ^P°,(r;R)  -J-  ^(r;  R)^ 


R)  = 


Hy(r;  R) 


1.15 


1.16a 


(A  /,R) 


^(r;  R) 


^(r;  R) 


r 


1.16b 


where  the  H^(r;  R)  are  the  adiabatic  state  wavefunctions  obtained  from  the  solution  of  eq.  1.5. 
The  g(J,  I,  R)  or  first  derivative  nonadiabatic  coupling  matrix  elements  which  (for  real-valued 
wavefunctions)  only  couple  PES's,  that  is  g(I,  I,  R)  =  0,  are  usually  the  quantities  of  principal 
concern  with  the  h(J,  I.  R)  and  k(J,  I,  R),  which  both  couple  and  modify  PES's  being  required  in 
full  quantum  mechanical  adiabatic  state  treatments. 
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(A)  Methodology 

Over  the  last  five  years  we  have  developed  a  unique  system  of  computer  algorithms39'40 
for  evaluating  g,  h  and  k  which  provide  significant  computational  advantages  over  previously 
existing41  methods.  We  employ  techniques  which  are  applicable  to  large  scale  Cl  wavefunctions 
(105  -10°  terms)  developed  from  state  averaged  MCSCF4^  reference  spaces  enabling  the  evaluation 
of  the  derivative  couplings  using  the  same  wavefunctions  which  yield  reliable  PES's.  This  is  an 
extremely  important  aspect  of  our  approach  since  the  nonadiabatic  coupling  matrix  elements  are  a 
sensitive  function  of  the  separation  of  the  PES’s  and  hence  of  the  nuclear  coordinates. 

We  derive  our  computational  efficiency  in  two  ways.  The  wavefunction  derivatives 
required  in  eqs.  1.15  and  1.16  are  evaluated  using  the  same  analytic  gradient  techniques  employed 
in  the  evaluation  of  energy  gradients.3913  Thus  we  obtain  savings  in  computational  effort  similar  to 
the  substantial  savings  achieved  when  analytic  energy  gradients  replace  energy  gradients  based  on 
divided  differences.  A  corollary  of  this  approach  is  a  unified  system  of  programs  for  obtaining 
both  derivative  couplings  and  energy  gradients.3913 

The  second  source  of  computational  efficiency  obtains  in  our  treatment  of  derivative 
couplings  attributable  to  overall  nuclear  translation  and  rotation.  For  Ra  an  overall  nuclear 
translation  or  rotation  the  derivative  couplings  matrix  elements  can  be  obtained  from  eqs.  1.15  and 
1.16  or  as  matrix  elements  of  electronic  operators  including  the  mass  polarization  operator  and  the 
total  electronic  angular  momentum  operator.40  By  replacing  the  nuclear  derivative  operator  in  eqs. 
1.15  and  1.16  with  these  (body  fixed  frame)  electronic  operators  the  costly  evaluation  of  orbital 
and  Cl  coefficient  derivatives  can  be  avoided.  Further  this  approach  permits  us  to  make  maximal 
use  of  molecular  point  group  symmetry  when  determining  the  adiabatic  states  an  advantage  which 
does  not  obtain  if  the  nuclear  derivative  formulation  is  used  exclusively.  Finally  the  observation 
that  for  overall  nuclear  translations  and  rotations,  alternative,  independent  computational  paths  to 
the  same  result  exist  -  one  exploiting  gradient  techniques  and  the  second  involving  only  matrix 
elements  of  electronic  operators  -  provides  a  stringent  test  of  the  precision  of  both  algorithms.  A 
detailed  discussion  of  these  points  can  be  found  in  reference  40. 
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(B)  Applications 

(i)  The  Energetic  System:  He  +  H2(Bl£^) 

In  the  following  we  consider  the  stability  of  a  proposed  high  energy  density  material  with 
respect  to  radiationless  decay.  This  study  provides  a  clear  demonstration  of  the  power  of 
theoretical  methods  to  (quickly)  assess  the  utility  of  a  proposed  energy  storage  system  by 
determining  the  topology  of  the  relevant  potential  energy  surfaces  (PES’s)  and  the  nonadiabatic 
interactions  which  couple  them. 

The  H2(B1Zy)  state  contains  over  11  eV  of  energy  relative  to  the  H2(X1X^)  ground 
state.32  It  has  been  suggested  that  the  charge  transfer  structure  (HeH)+  -  H“  which  consists  of 
two  stable  moieties,  HeH+  and  H~  might  exist  as  a  stable  species  on  either  the  2^'  PES  (as  a 
maximum  ionicity  excited  state  -  MIES43  structure)  or  on  the  ground  1 1  A'  PES  and  thereby 
provide  a  chemical  means  of  storing  the  electronic  energy  in  the  B*X^  state  of  H2. 

To  address  this  question  the  1,  2^'  states  of  the  He-H2  system  were  characterized  using 
ab  initio  MCSCF/second  order  Cl  wavefunctions.44  These  surfaces  correlate  with  the  He('S)  + 
H2(X1Zg,  B1!^)  system  states.  This  work  provided  the  first  extended  treatment  of  the 
nonadiabatic  interactions  between  the  1, 2*A'  states  reporting  first  derivative  nonadiabatic 
couplings  (eq.  1.14)  at  over  100  nuclear  configurations.  Two  extrema  on  the  2*A’  PES  were 
determined  at  the  Cl  level,  a  global  minimum  and  a  saddle  point.  In  each  case  the  character  of  the 
extremum  was  determined  from  the  eigenvalues  of  the  second  derivative  matrix.  While  this  work 
built  on  previously  reported  theoretical  treatments  of  these  PES's,45  both  qualitative  and 
quantitative  differences  with  previous  work  were  found.  The  predicted,  entrance  channel,  saddle 
point  on  the  2lA'  PES  gave  a  barrier  of  1.5  Kcal/mol  which  is  significantly  lower  than  previous 
estimates.453  The  global  minimum  on  the  2^'  PES  corresponds  to  a  charge  transfer  structure, 
(HeH)+-H~  and  is  stable  by  38.99  Kcal/mol  relative  to  the  He(!S)  +  ^(B1!^)  asymptote.  The 
structure  and  energy  of  this  extremum  are  in  qualitative  accord  with  the  previous  ab  initio  results  of 
Farantos  et  al.43 

The  most  significant  new  finding  in  this  work  was  the  demonstration  of  the  existence  of  a 
'seam'  for  which  the  1  lA'  and  2*A'  PES's  are  nearly  degenerate,  AE  =  E^A')  -  E^A')  < 


14 


0.001-0.5  Kcal/mol.  This  extended  region  of  large  nonadiabatic  effects  occurs  for  general  Cs, 
rather  that  C2v  or  G»v  geometries,  and  includes  the  region  of  the  minimum  on  the  2*A’  PES  noted 
above.  This  feature  of  the  1 ,  2 1  A’  PES's  should  significantly  influence  the  quenching  of 
H2(B1X^)  by  He  and  alter  expectations  of  the  lifetime  of  the  state  corresponding  to  the  global 
minimum  on  the  2*A'  PES.46  It  is  anticipated  that  the  lower  entrance  channel  barrier  noted  above 
and  the  near  degeneracy  seam  reported  here  will  serve  to  increase  the  calculated  quenching  cross 
section47  which  is  appreciably  lower  than  the  measured  value.4**  The  large  nonadiabatic  coupling 
between  the  1 ,21A’  wavefunctions  in  the  vicinity  of  the  global  minimum  on  the  2*A’  PES  will 
serve  to  decrease  the  predicted  lifetime  of  this  'state'.  To  facilitate  dynamical  studies  of  these 
questions  first  derivative  nonadiabatic  coupling  matrix  elements  were  determined  in  the  vicinity  of 
the  'seam'. 

The  possible  stability  of  the  charge  transfer  structure  (HeH)+-H~  on  the  1  *A'  PES  was 
considered  using  a  surface  walking  technique.  In  the  vicinity  of  the  above  noted  seam  an  exchange 
of  the  charge  transfer  structure  between  the  1, 2JA'  states  is  observed.  Surface  walks  on  the  1  !A’ 
PES  which  follow  the  energy  gradient  from  nuclear  configurations  corresponding  to  the  charge 
transfer  structure,  (HeH)+-H~,  were  performed.  In  no  case  was  a  stable  triatomic  structure 
corresponding  to  a  charge  transfer  configuration  found. 

C.  Avoided  and  Allowed  Surface  Crossings 

(i)  Efficient  Determination  of  Regions  of  (Avoided)  Surface  Crossings  Using  Analytic  Energy 
Gradients 

A  key  to  the  understanding  of  electronically  nonadiabatic  processes  is  the  determination  of 
the  regions  of  nuclear  configuration  space  for  which  the  potential  energy  surfaces  in  question  are 
closely  spaced  or  actually  cross.  It  is  in  these  regions  that  nonadiabatic  effects  are  most  important. 
The  location  of  such  regions  of  coordinate  space  can  be  quite  time  consuming.  Recent  advances  in 
analytic  gradient  techniques49"53  have  reduced  considerably  the  labor  required  for  determining 
single  surface  extrema.  As  part  of  this  research  effort  we  have  begun  the  development  of 
algorithms  to  bring  comparable  efficiencies  to  the  determination  of  (avoided)  surface  crossings.54 
These  algorithms  utilize  the  same  state  averaged  MCSCF/direct  symbolic  matrix  Cl  wavefunctions 
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which  are  used  to  characterize  the  spin-orbit  and  spin-spin  interactions.  They  are  based  on  analytic 
gradient  techniques  thereby  achieving  the  considerable  efficiencies  already  available  to  us  in  the 
location  of  extrema  on  a  single  potential  energy  surface. 

This  project  was  motivated  by  our  previous  studies  of  spin-forbidden  and  electronically 
nonadiabatic  decay  pathways  in  polyatomic  energetic  materials  As  noted  above  we  have 
considered  (1)  the  mechanism  of  the  spin-forbidden  radiationless  decay  of  hydrazoic  acid 
N3H(X1A')  — >  N2  +  NH(X3L*)  and  (2)  the  formation  of  the  [unstable  with  respect  to 
decomposition  into  Ke+H2(1S^)]  charge  transfer  state  HeH+-H-  from  He  +  H2(BJ£*).  In  order  to 

o 

characterize  the  electronic  structure  aspects  of  process  (1)  it  was  necessary  to  determine,  in  addition 
to  the  requisite  spin-orbit  coupling  matrix  elements,  the  minimum  energy  crossing  point  of  the 
1 1A'  and  l^A"  potential  energy  surfaces.  Similarly  the  elucidation  of  the  mechanism  of  process 
(2)  required  the  determination  of  the  nonadiabatic  coupling  matrix  elements  between  the  1  'A'  and 
2*A'  potential  energy  surfaces  as  well  as  a  seam  of  avoided  crossings  [R(r),  y(r)]  where  r  is  the 
H2  distance,  R  is  the  He-H2  distance  and  y  is  the  He-H2  angle.  In  the  above  noted  studies  the 
minimum  energy  crossing  point  and  seam  of  avoided  crossings  were  determined,  somewhat 
tediously,  by  consideration  of  the  individual  potential  energy  surfaces.  The  analytic  gradient 
driven  algorithms  for  determining  allowed  and  avoided  crossings  of  potential  energy  surfaces 
would  have  reduced  considerably  the  computational  effort  required  to  treat  these  problems. 

These  procedures  have  broad  applicability.  In  particular  these  algorithms  can  be  expected 
to  find  considerable  use  in  the  characterization  of  high  energy  density  materials  as  these  species  are 
frequently  characterized  by  locally  stable  regions  of  a  potential  energy  surface  which  are  isolated 
from  more  stable  structures  by  (avoided)  surface  crossings. 
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